rm(list=ls())


load('directory of Res_table2')


### Minneapolis
ce.min=apply(ce.min.rho0[,1:5],2,mean)   

ce.min.ci=round(apply(ce.min.rho0[,1:5],2,quantile,probs=c(0.025,0.975)),3)


##### Denver
ce.den=apply(ce.den.rho0[,1:5],2,mean)   

ce.den.ci=round(apply(ce.den.rho0[,1:5],2,quantile,probs=c(0.025,0.975)),3)



#### Pooled
p=562/955

ce.pool=ce.den*p+ce.min*(1-p)
ce.pool.ci=  p*ce.den.ci[,]+ (1-p)*ce.min.ci[,]



